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ABSTRACT 

Context. High energy electrons accelerated during solar flare are abundant in the solar corona and in the interplanetary space. 
Commonly, the number and the energy of non-thermal electrons at the Sun is estimated using hard X-ray (HXR) spectral obser- 
vations (e.g. RHESSI) and a single-particle collisional approximation. 

Aims. To investigate the role of the spectrally evolving Langmuir turbulence on the population of energetic electrons in the solar 
corona. 

Methods. We numerically simulate the relaxation of a power-law non-thermal electron population in a collisional inhomogeneous 
plasma including wave-particle, and wave-wave interactions. 

Results. The numerical simulations show that the long-time evolution of electron population above 20 keV deviates substantially 
from the collisional approximation when wave-particle interactions in non-uniform plasma are taken into account. The evolution of 
Langmuir wave spectrum towards smaller wavenumbers, due to large-scale density fluctuations and wave-wave interactions, leads to 
an effective acceleration of electrons. Furthermore, the time-integrated spectrum of non-thermal electrons, which is normally observed 
with HXR above 20 keV, is noticeably increased due to acceleration of non-thermal electrons by Langmuir waves. 
Conclusions. The results show that the observed HXR spectrum, when interpreted in terms of collisional relaxation, can lead to an 
overestimated number and energy of energetic electrons accelerated in the corona. 
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1. Introduction 

Solar flares provide a number of theoretical challenges for var- 
ious aspects of electron dynamics, including their accelera- 
tion and transport. One of the important aspect of solar flares 
is the high acceleration efliciency of electrons, seen via their 
Hard X-ray emission (HXR) in the dense atmosphere foot- 
points (e.g.^rown et al. 1990; Shibata 1999; Aschwanden 2002^ 
pn et al. 2003; Lin 2006; Bian et al. 2010; Battaglia et a l.ll201U 
Battaglia & Kontar. ,2011: .Christe et al.. ,2011: .Krucker et alJ 



transport. In this approximation the energy of an electron E - 
mv^/2 decreases with time due to binary Coulomb collisions, so 
that for non-thermal particles v(f) » Vtc, one can write 



dv 
dt 



(1) 



201 II) . It was realized early on th at the solar corona can store 
large amounts of magnetic energy (ISweetlll958l;lPriest & ForbesI 
20021) . that can be released to increase the kinetic energy of the 



• • i surrounding plasma. While the exact location and properties of 
. ^ ' electron acceleration site are still under debate, spatially re- 
k> solved HXR observations, notably with RHESSI, point to the 
5_] coronal origin of electron acceleration and subsequent trans- 
port of energetic electrons towards the d enser lower corona (e.g. 
' iHolman et all 1201 ll: iKontar et al.ll201 ll as the recent reviews). 
HXR emission from the dense regions is then used as the model 
dependent diagnostic of electron acceleration in the corona. The 
low dynamic range of current hard X-ray instruments is insuf- 
ficient to observe the evolution of the energetic electron distri- 
bution enroute between the coronal and the footpoint parts. The 
exception to this rule is probably the radio observations of emis- 
sions from energetic electrons (e.g. Fleishman et al. 2011). 

The propagation of energetic electrons between the coronal 
and the footpoint sources is often treated as a simple collisional 



where F = 47re'*n In A/m^, n is the number density of the sur- 
rounding plasma, e is the electron charge, m is electron mass, Mje 
is the electron thermal velocity, and In A ^ 20 is the Coulomb 
logarithm. This treatment obviously ignores all collective inter- 
actions in plasma, but gives a simple relation between the ac- 
celerated and the resulting mea n electron spectra re sponsible for 
the observed HXR spectra (e.g. lBrown et al. 2003lPl 

Realizing the importance of various streaming instabilities 
and the subsequent wave-particle interactions, a number of mod- 
els have included the generation of Langmuir waves. Indeed, 
Langmuir wave turbulence can be effecti vely generated when 
the bump on tai l distribution is formed ( iKaplan & TsvtovichI 
I1973L iTsvtovichI [1995 ). Firsdy, the initiaUy stable 5//5v < 
electron distribution function /(v) can be turned unsta- 
ble i.e. dfldv > 0, when the fast electrons overtake the 
slower ones. In the case of collisional relaxation of the elec- 
tron population, since slower particles lose their energy faster, 
Tcoi{E) ^ mv-'/(2r), a "gap" distribution can appear from 
an initially power-law distribution on collisional time-scales. 
This "gap" distribution of electrons creates a positive slope 
df/d v > 0, ( Svi'ovatskii & Shmeleva 1972; Emslie & Smit^ 
il984l:lHamilton & Petrosianll 19871) . which could lead to growth 
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of Langmuir waves. The beam-plasma instability acts to de- 
crease the gradient of the electron distribution in velocity space, 
while Coulomb collisions tend to restore it. The net result is 
a plateau-like distribution, with the plateau height slowly de- 
creasing as energy is lost to the background plasma v i a col- 
lisions. It has been shown (e.g. [Ham ilton & PetrosianI 119871; 
lMcClementslll987t iHannah & Kontaii 1201 1 ) that even for the 
time independent injection of electrons into plasma, the elec- 
trons should effectively generate Langmuir waves, however the 
time or space integrated electron spectrum (to compare with 
the HXR observations) is only weakly affected by this pro- 
cess. The overall evolution of the beam-plasma system is com- 
plicated by various non-linear proces ses, which can a ffect the 
spectrum of Langmuir turbulence (see lTsvtovichlll995L as a re- 
view). Due to these processes, the energy of Langmuir waves 
can be transferred from one range of phase velocities to another, 
hence suppressing the growth of Langmuir waves. The removal 
of Langmuir waves out of resonance with the electron beam can 
potentially lead to nonlinear stabilization of the beam-plasma 
instability ("Papadopoulos' 1975; Vlahos& Papadopoulosl l 1 9791 
[Rowland & Vlahos 1985; Muschietti & Dum 1991)'. While these 
earlier works were predominantly analy tical estimates, the 
recent weakly turbulent treatments (e.g. IZiebell et al.] 120011; 
iKontar & Pecsell 120021; Iziebell et all 1201 lb are self-consistent 
and capture a complex interplay between wave-particle and 
wave-wave processes in collisionless plasma. For beam relax- 
ation in a plasma, these processes can lead to the appearance 
of electrons with velocities above the maximum injected. At 
shor t scales, beam-plasma interaction Vlasov or PIC simulations 
(e^.'Kariickv & Kaspa roval2009l;lTsiklauril2010l;lDaldorff et alj 
1201 1 ; Tsiklauri 2011) are particularly useful to capture phase- 
space effects of beam-plasma interaction at short time scales. 
However, the interpretation of solar flare HXR emission requires 
knowledge of electron evolution at much longer time-scales, in- 
accessible by Vlasov or PIC approaches. 

In this paper, we investigate the evolution of non-thermal 
electrons and Langmuir waves using weak-turbulence theory in a 
coUisional, inhomogeneous solar coronal plasma. We show that 
the evolution of the Langmuir wave spectrum caused by plasma 
inhomogeneities and wave-wave processes leads to effective ac- 
celeration of high energy electrons. The resulting HXR spectrum 
is found to be strongly affected by the evolution of Langmuir tur- 
bulence, so that if the HXR spectrum is interpreted in terms of 
a collisional model, the accelerated electron spectrum will be 
overestimated. 



2. Weak turbulent evolution of energetic electrons 
in a collisional plasma 

The weakly turbulent collisional relaxation of an energetic elec- 
tron population is considered in a collisional plasma typical for 
the solar corona. The evolution of energetic electrons is de- 
scribed using weak turbulence theory including wave-particle 
and wave-wave inter actions in non-uniform plasma. The corre- 
sponding equations ( Vedenov & Velikho"vlll963l ; IVedenov et alj 
11 9671; iTsvtoviciill 19951) ' governing the electron distribution func- 
tion /(v, /) [electrons cm"^ (cm/s)"' ] and the spectral energy 
densities of Langmuir Wk, and ion - sound waves W'^ [ergs cm"^ 
] are as in lKontar & Pecselil (l2002h : 



of ov \ V ov / 



dt 



dajpeix) dWk 



-'f 

dx dk 

vln 



nk^ ov 



/ + StdecayiWk, W'^. 



(3) 



Here, Sideca\(Wk, W^) is the integral describing the parametric 
decay of a Langmuir wave with an ion-sound wave represented 
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(5) 



where the spectral energy density is normaUzed such, that 
J Wk dk is the energy density of the waves [erg cm"^]. The sys- 
tem of equations (|2]|5]l describes the weakly turbulent evolution 
of electrons and Langmuir waves in the presence of density inho- 
mogeneities and/or ion-sound waves. The spontaneous emissio n 
of Langmuir wav es has been treated as in iHannah et al ] d2009h ; 
iReid et al] (1201 ll) and the sys tem was numer ically integrated us- 
ing a numerical scheme as in iKontaJ (1200 1 d) . 

As we are interested in the evolution of the system 
at time-scales t » T cq/, the collisional operator (e.g. 
iLifshitz & Pitaevskiilll981h is included to account for the binary 
collisions 



pi 

si.jf) = r— 



'4. 



'Te 



dv] 



(6) 



where the first RHS term describes the systematic drag on 
energetic particles and the second term diffusion in velocity 
space due to binary collisions. We note that the collisional drag 
is larger than energy losses due to spontaneous emission of 
Langmuir waves [second term in the right hand side of Eq. ©I 
by a non-thermal electron in a plasma. In other words, only 
a small part of the energy lost by an electron with velocity v 
goes into spontaneously generated Langmuir waves given by Eq. 
Q. Finally, the collisional damping rate for Langmuir waves 
Jcoi - r/4v^^ is added to Eq. ([3]l. 

2. 1. Initial electron distribution function 

Let us consider the evolution of a non-thermal electron popula- 
tion with the initial electron distribution go{v) in a Maxwellian 
plasma (Fig.[TJ: 



/(v,f = 0) = 



V2^' 



■ exp 



,2 \ 



+ ^o(v), 



(7) 



where the non-thermal particle distribution ^o(v) is initially a 
power law /(v, f = 0) ~ v"^'* for v > v^ = lOvje and flattens at 
low velocities v < IOvtc', 



(2) §o(v) 



2nb r(6) 
V^Vfe r{6- i) 



[l+(v/v,)2]" 



(8) 
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where 6 is the power law index for the energetic particles in 
energy space, nh the number density of non-thermal electrons, 
nt, « n, and r(x) denotes the gamma function. The initial elec- 
tron distribution is normalized to the electron number density 
[electrons cm""*], so that 



Time, s 



f 

Jo 



^o(v)dv = nt. 



(9) 



The initial level of Langmuir waves can be calculated assuming 
an equilibrium with Maxwellian electron distribution function, 
ignoring all non-linear terms and setting dWk/dt = in Eq. (O 



W(k, t^O) 



1+ A/i^^^^4.exp( 



khT, 



47r2 



1 + 

1 67m. 



Vf^'exp( 



(10) 



where kh is Boltzmann constant, Tp = mWj^ the electron temper- 
ature of the plasma, and Aoe - ^Te/'^pe is the Debye length. The 
equation ([TOl) can be reduced to W(k,t - 0) - j^k^ln(jj^) 
in the coUisionless limit for ycoi — > 0, which is the thermal 
level of plasma waves in a coUisionless M axwellian plasma 
dKaplan & TsvtovichllT97l lTsvtovichllT995h . For the problem 
considered, the exact initial level of plasma waves is not im- 
portant, as the governing equations quickly establish a balanced 
level, which is a few orders of magnitude less than the level 
driven by the instability. 



3. Evolution of the electron distribution function 

We assume a background plasma similar to the solar corona, 
with plasma frequency of fj,g - Upelln - 2 GHz and an elec- 
tron and ion temperature of T,- = = 1 MK. The correspond- 
ing plasma density is then 5.0 x 10 '"cm"^ and the coUisional 
timescale is TcoI = I/Tco/ = 7 x 10"^ seconds. For the electron 
beam we take rih/n = 1 x 10"^, and (5 = 4, giving a power-law 
index of 8 in velocity space above approximately v;,. 

3.1. CoUisional relaxation 

Under the influence of Coulomb collisions, the electron distribu- 
tion function for v » Vxe evolves according to 



dt dv\dt^j dv\ v2 



(11) 



which can be straightforwardly solved for arbitrary initial dis- 
tribution function, goiv) so that f{v,t) = [\/u(v,t)]^go{u(v,t)), 
where m(v, f) = (v-' + SFf)'^"'. The time integrated, or mean elec- 
tron distribution, will be the so-called thick-target electron spec- 
trum. This expression is, as expected, sim ilar to the space inte- 
grated coUisional thick-target model (e.g. iDubovl [19631 iBrownl 
[197 1; Brown et al. 2003) and directly corresponds to the elec- 
tron spectrum responsible for X-ray flux. Figure ([TJ shows the 
time evolution of the electron distribution and the time inte- 
grated electron flux spectrum. The time-dependent positive slope 
df/dv > is formed for velocities v < (3rO'^^ while the tail at 
V > (3Ff)'^"^ is weakly affected by collisions. As time progresses, 
a larger range of electron velocities shows a distribution function 
which grows with velocity. 
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Fig. 1. CoUisional relaxation of an electron beam in a plasma. 
Top and middle panels; the electron distribution function /(v) 
and the spectral energy density of Langmuir waves (here fixed at 
the thermal level), respectively. Each coloured line shows the 
distribution at a different time. Bottom panel; the time aver- 
aged electron flux spectrum [electrons keV"' cm"^ s"'] plotted 
against electron energy. The dashed line shows the flux if the 
beam remained the initial power-law. 



3.2. Langmuir wave generation in uniform plasma 

In the case of a uniform plasma (dcope/dx = 0) and for sim- 
plicity ignoring all wave-wave terms, the coUisionally formed 
unstable distribution will lead to instability and effective gen- 
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eration of plasma waves (Fig. |2]|. The generation of Langmuir 
waves quickly flattens the distribution function and forms a 
plateau between ~ 'i'Vje < v < (SFf)'^"' (Fig. |2|. Over time, 
more electrons are lost due to collisions, and the plateau in 
velocity space becomes wider. As soon as the beam density 
has dropped, Langmuir waves are no longer generated appre- 
ciably, and the evolution of the electron beam becomes again 
purely collisional. The important aspect is that while the instan- 
taneous distribution of electrons in the cases of quasi-linear and 
collisions-only relaxation are very different, the time-integrated 
spectrum of energetic electrons is almost identical. This con- 
firms the pre viously published results (Hamilton & Petrosianl 
1987HMcCle"ments, 1987: Hannah et al, 2009; Hannah & Kont3 



Time, s 



20111) that the quasi-linear interaction weakly affects the time- 
integrated spectrum. The energy is quickly transferred from elec- 
trons to the plasma waves and back, the time integrated amount 
of energy for a given r esonant velocity is almost constant (e.g. 
IVedenov & Velikhovl 1 963 ■ MeFNik et al.l 1 999) . The colHsional 
energy loss of Langmuir waves does not faciUtate the transfer of 
energy between different velocities, but simply reduces the en- 
ergy. 



3.3. Spectral evolution of Langmuir turbulence 

The Langmuir wave spectrum can quickly evolve in k space, 
changing the wave energy available at different resonant veloc- 
ities, V - cope/k. Firstly, the nonlinear wave-wa ve interactions 
change the distribution of Langmuir wave energy (Vedenov et al.l 
1967; PapadoDoulos et al. 1974; Goldman & Dubois 1981 
Barta & Karlickv 2000; Y oon etal.l 120051: jFouquet & Pesmg 
2008; Palastro et al. .2009.) and secondly, the waves can be 
effectively scattered or refracted by inhomogeneous plasma 
(Rvutov '1969"; 'Breizman & Ruvtov '1969; 'Krasovskii' '19781 
Coste et al. 1975; Nishikawa & Rvutov 1976; Smith & Simd 
1979t j Asadulfin et al. 1980; Mekose & Cramei] [l989t FKonteJ 
2001al[briLi et al..,2006; ,Ziebell et al..,2011l) . 



3.3.1 . Inhomogeneous plasma - constant density gradient 

To illustrate the role of the spectral evolution of Langmuir waves 
on the electron distribution function, we consider the relaxation 
of energetic electrons in a plasma with a constant density gra- 
dient, so that the term describing refraction of Langmuir waves 
is 



dojpe 

dx 
where 



L 



dx 



(x) I dnpeix) 



dx 



(12) 



(13) 



is the characteristic scale of density inhomogeneity. Langmuir 
waves propagating into the region of higher/lower density expe- 
rience change in wavenumber lo k + Ak, with Ak ^ +(L)peAt/\L\, 
depending on the sign of the density gradient. This in turn affects 
the resonant velocity v = cuipe/k, so electrons of smaller/larger 
velocities can now interact with this wave. For electrons prop- 
agating into denser regions, the plasma density grows, causing 
generated Langmuir waves to have a negative shift in k, Ak ^ 
-CL)peAt/\L\. Therefore, Langmuir waves generated at larger k are 
shifted towards smaller k (higher phase velocities), and hence, 
can be effectively absorbed by faster electrons with df/dv < 0, 
which are consequently accelerated. For a density gradient of 
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Fig. 2. The same as in Figure [T]but for collisional relaxation in- 
cluding Langmuir wave generation in homogeneous plasma. In 
the bottom panel, the black line shows the collisions only case, 
the blue included wave generation. 



L > 10^ cm, the waves are shifted faster than the characteristic 
time of collisional relaxation, so for L ~ ^TeJ^-li ~ 10^ cm, sub- 
stantial energy can be transferred to higher electron velocities. 
Figure [3] shows clear acceleration of energetic electrons above 
~ 3QvTe as the Langmuir waves shift to higher phase velocities. 
Importantly, these accelerated electrons are evident in the time 
averaged spectrum, hence the resulting Hard X-ray spectrum 
will be affected and the role of Langmuir waves is detectable. 
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Fig. 3. The same as in Figure |2] but for a plasma with constant 
density gradient as in Equation fT2l with L - 10* cm. 



3.3.2. Density fluctuations 

However, even in the absence of a constant density gradient, 
the Langmuir wave spectrum is also expected to evolve. For in- 
stance, the solar coronal plasma, as any real plasma, has a fluc- 
tuating part, n + 6n{x, t). For large scale density perturbations 
1 /L << (Dpelvb, the resonant conditions for Langmuir wave de- 
cay or scattering cannot be satisfied and the Langmuir waves 
will be scatter ed with a sniall ch ange in k (Nishikawa & Rvutov 
[1926; Goldm an & DuboisllTgsl . Let us consider density fluctu- 
ations with {(6n/n)^) « 1, zero mean (Sn) = and a Gaussian 



spectrum |5n^|n,^ ~ sw(~1^/lo) expC-Q^/Qg), where qo and Qo 
are their characteristic wavenumber and frequency. Under the in- 
fluence of these turbulent density fluctuations, the spectral evo- 
lution of Langmuir waves in k becomes diffusive 

where () denotes ensemble average over density fluctuations. 
The diffusion coefficient describing this random refraction of 
waves is 

D,^fulAMnfff^[l.'^f'\ (15) 

where = 2uoojpe/3Vj^ and mq = Oq/^o- Assuming a modest 
level of turbulence 6n/n - 10^^, with a characteristic speed close 
to the speed of sound in the corona uq = fio/^o ~ 10^ cm s ' and 
long wavelength l/qo ~ 10"* cm, this spetral diffusion D^, sub- 
stantially changes the spectrum of Langmuir waves in a few col- 
lisional times. Figure |4] shows the evolution of the electron dis- 
tribution function and the spectral energy density of Langmuir 
waves. Langmuir waves diffuse in A:-space, from an initial value 
of A: ~ ojpe/vh, so in time the spectrum of Langmuir waves be- 
comes broader as time increases. The absorption of Langmuir 
waves by electrons above lOvje leads to electron acceleration. 
As in the constant gradient case, the integrated electron flux 
for the fluctuating density case changes considerably at energies 
above ~ 20 keV. 

3.3.3. Tliree-wave interaction of Langmuir waves 

For density fluctuations with wavenumber q ~ copelvb the dif- 
fusive treatment is no longer valid, and we must individually 
consider the interactions between Langmuir waves and other 
modes in the plasma. In particular, we are interested in those 
modes which may be excited by the high level of Langmuir 
waves present due to the beam-plasma interaction. Following the 
treatments in uniform dZ iebell et al. 2001) and inhomogeneous 
(iKontar & Pecselil 120021 : iZiebell et al.. i201 L) plasmas, we now 
consider the three wave process, I'; I ^ I' + s, where / and /' are 
Langmuir waves and s is a sound wave. This process efficiently 
scatters electron beam generated waves with wavenumber k, into 
secondary Langmuir waves with wavenumber k' x -k. The 
decay of a Langmuir wave generates a ion-sound wave with 
q ^ 2k. Every scattering decreases the absolute value of the 
Langmuir wave number hy k*^ - 2 V'«e/'«i Vl~+3777TJ/(3/loe). 
Hence, repeated wave-wave processes cause Langmuir waves to 
appear at lower k (larger phase velocities). As in the previous 
cases, Langmuir waves can be absorbed by the tail of the distri- 
bution, leading to acceleration of non-thermal electrons. 

The simultaneous evolution of both io n- sound waves and 
Langmuir waves has been already simulated (IZiebell et al.l200ft 
IKontar & Pecselil 120021: Iziebell et al.ll20Tll e.g.) and does pro- 
duce the effect outlined in the previous paragraph. Here we limit 
ourself to a more simple but quite illuminating situation with a 
fixed level of density fluctuations produced by ion-sound waves. 
We assume the level of ion-sound waves to be thermal. 

The results are presented in Fig. |5] As in the cases of con- 
stant density gradient and random density fluctuations, we see 
a time evolving plateau distribution up to ~ 6QvTe at f - Is. 
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Fig. 4. The same as in Figure [2] but for the case with random Fig. 5. The same as Figure |2] but with ion-sound wave interac- 
density fluctuations, given by Equation ( fTSl l. tions as well as density fluctuations. 



The nonlinear decay term has a strong impact on the distribu- 
tion above SOvje due to the high level of waves for k < koe/^O. 
As the wave-wave terms are particularly effective in shifting the 
Langmuir wave energy to low k, making Langmuir waves more 
isotropic, the number of accelerated electrons at energies above 
100 keV also increases. Below 100 keV, the effect in the time 
integrated spectrum is similar to the previous cases. 



4. Discussion and conclusions 

The above results show the significant effect that Langmuir 
waves have on the corresponding electron distribution function. 
A constant density gradient, or random long-wavelength density 
fluctuations, produce enhancement of the electron distribution in 
the range 30-60 Wje- On the other hand, the wave-wave interac- 
tions with short wavelength fluctuations q ~ 2k show a more pro- 
nounced effect at velocities v ^ 60vTe- Unlike the case of a uni- 
form plasma, the time averaged electron spectrum is noticeably 
different from collisional relaxation. As the time averaged spec- 
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-1000 



Energy, keV 

Fig. 6. The effective energy loss rate of an electron {dE/dt)eff 
using Equation (fTTT i. at 0.5s. Collisional losses {dEldt)eff - 

—T yjn^JlE^^ from Equation[T](black dashed line), collisional 
losses as in Figure[T](black line); Constant density gradient as in 
Figure [3] (blue line); Random density fluctuations as in Figure |4] 
(red line); Non-linear interactions in inhomogeneous plasma as 
in Figure |5] (green line). 

trum can be directly compared to Hard X-ray observations, the 
role of Langmuir waves could be observable. The resulting HXR 
spectrum above ~ 20 keV is higher than expected from colli- 
sional relaxation, and if interpreted as thick-target relaxation, the 
number of energetic electrons injected is overestimated. 

Given the evolution of the electron distribution, /(v, f), we 
can calculate the effective energy change rate. Evidently the ac- 
tual energy loss/gain rate is given by the kinetic equations (l2]|5]l, 
but assuming a single particle evolution often used in the lit- 
erature (e.g. Brown et al. 2009) the effective energy loss rate 
{dE/dt}eff can be defined as the one satisfying a continuity 
equation: 







(16) 



where f(E, t) - mw(E)f(y(E), t) is the energy distribution func- 
tion. The effective energy loss/gain rate satisfying this equation 



beams (e.g. lHannah & Kontarll20lTl) . when electrons are accel- 
erated in bunches rather as a continuous stream as expected in 
some acceleration models. The acceleration of energetic elec- 
trons in non-uniform plasma is due to Langmuir waves gener- 
ated by the non-thermal electrons themselves and does not in- 
crease the total energy of the beam-plasma system. On the con- 
trary, the energetic electrons lose their energy not only via colli- 
sions but also via generation of Langmuir waves and their sub- 
sequent absorbtion via collisions. Therefore, although the low 
energy electrons lose their energy faster, the higher energy elec- 
trons are gaining energy and hence can generate HXR emission 
for a longer time. This "self-acceleration" of some electron pop- 
ulation reduces the number of accelerated electrons necessary to 
explain the observed HXR fluxes. This process also produces a 
similar effect to a re- acceleration from ex ternally applied electric 
fields as discussed bv lBrown et al.l (l2009l) . but evidently does not 
require additional energy input. 

A few important comments are due on the application of 
these results to the interpretation of solar flare X-ray spec- 
trum. The overall effect of electron re-acceleration due to evo- 
lution of Langmuir turbulence is beam-density dependent. The 
higher the beam density, the stronger the generation of Langmuir 
waves, and therefore the effects considered should be more 
evident in HXR data. Therefore, the stronger the flare, the 
larger the overestimation of non-thermal electrons, when inter- 
preted by a collisional thick-target model. Further, these sim- 
ulations need to be extended to account for the spatial trans- 
port of energetic electrons in solar flare conditions, a relativis- 
tic tr eatment of the electron s included, and the role of whistlers 
(e.g. Steoanov & Tsap 2002) and kinetic Alfven waves (e.g. 
Bian et al. 2010) considered. The relaxation produces an en- 
hanced level of La ngmuir waves, which could lead to plasma mi- 
crowave emission (Ems lie & SmithI 19841; [Hamilton & PetrosianI 
1987) and needs to be investigated further. However, as the spec- 
trum of Langmuir turbulence evolves with time, the wave-wave 
and density gradient effects need to be included. 

In summary, the evolving Langmuir turbulence spectrum has 
a noticeable effect on both the instantaneous and time-averaged 
distribution functions leading to the overestimation of electron 
flux (and hence energy) in non-thermal electrons > 20 keV, if 
interpreted in terms of a standard collisional thick-target model. 
Therefore, if Langmuir waves are effectively generated in flares 
and their spectrum evolves, the number of accelerated electrons 
could be substantially smaller than inferred from HXR thick- 
target interpretation. 



1 



f{E,t) 



r 



df{E',t) 
dt 



dE'. 



(17) 



If the evolution of particles is governed by collisional only, the 
effective loss rate is simply {dEldt)eff - -T ^fn^JilE) (Fig. 
|6ll. The evolving Langmuir wave turbulence notably changes 
{dE ldt)i,ff, decreasing the effective energy loss at energies of 
a few keV, increasing the loss rate in the intermediate range near 
10 keV, and leading to acceleration and/or reduced energy loss 
rate at tens of keV and above. This effective energy change rate 
increases the flux of electrons above 20 keV by a factor of a few 
to a few tens (shown in Figures |3]|5]l. 

We emphasize that the acceleration of energetic electrons 
is due to spectral evolution of Langmuir waves. The Langmuir 
wave generation in a uniform plasma only weakly changes the 
time integrated HXR spectrum and, hence, is not likely to be 
detected in current HXR measurements. The important excep- 
tion is the generation of waves by spatially localized electron 
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